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The Green-function approach has been extended to the moderate baryon 
density region in the framework of an extended Nambu-Jona-Lasinio model, 
and the thermodynamic potential with both chiral and diquark condensates 
has been evaluated by using the massive quark propagator. The phase struc¬ 
ture along the chemical potential direction has been investigated and the 
strong competition between the chiral and diquark condensate has been an¬ 
alyzed by investigating the influence of the diquark condensate on the sharp 
Fermi surface. The influence of the diquark condensate on the quark prop¬ 
erties has been investigated, even though the quarks in the color breaking 
phase are very different from that in the chiral breaking phase, the difference 
between quarks in different colors is very small. 
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I. INTRODUCTION 


QCD phase transitions along the baryon density direction has attracted much attention 
recently since it was found that the color-superconducting gap can be of the order of lOOMeV 
and 1^, which is two orders larger than early perturbative estimations [Q. However, till 
now, there is no uniform framework discussing the phase structure in the wide region of 
chemical potential fi from about SOOMeU to lO^MeV. 

In the idealized case at asymptotically high baryon densities, the color superconductivity 
with two massless flavors and the color-flavor-locking (CFL) phase with three degenerate 
massless quarks have been widely discussed from hrst principle QCD calculations, see @] and 
references therein. Usually, the diagrammatic methods are used in the asymptotic densities. 
The Green-function of the eight-component held and the gap equation were discussed in 
details in i i- Neither current quark mass nor chiral condensate are necessary to be 
considered because they can be neglected compared with the very high Fermi surface. At 
less-than-asymptotic densities, the corrections of non-zero quark mass to the pure CFL phase 
can be treated perturbatively by expanding the current quark mass around the chiral limit 


i i- 

For physical applications, we are more interested in the moderate baryon density region, 
which may be related to the neutron stars and, in very optimistic cases even to heavy-ion 
collisions. Usually, effective models such as instanton, as well as Nambu-Jona-Lasinio (NJL) 
model, are used. The model parameters are hxed in the QCD vacuum. In this region, 
the usual way is to use the variational methods working out the gap equations from the 
thermodynamic potential [0- except for and 13 in the instanton model, where 
the quark propagator was evaluated explicitly, but the form is complicated. 

One of our main aims in this paper is to apply the Green-function approach in the 
moderate baryon density region. To work out the phase structure from hadron phase to 
the color superconducting phase, one should deal with the chiral condensate and diquark 
condensate simultaneously. Because the chiral condensation contributes a dynamic quark 
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mass, it is not reasonable any more in this density region to treat the quark mass term 
perturbatively, like in |^. By using the energy projectors for massive quark, we will evaluate 
the Nambu-Gorkov massive quark propagator, which will help us deal with the chiral and 
diquark condensate simultaneously. 

In the normal phase, the quarks in different colors are degenerate, while in the color 
breaking phase, it is natural to assume that the quarks involved in the diquark condensate 
are different from that not participating in the diquark condensate. In and [l^, different 
masses for the quarks which participate and not participate in the diquark condensate are 
introduced. However, we will see it difficult to get the mass expression for the quarks 
participating in the diquark condensate, because the particles and holes mix with each other 
and the elementary excitations are quasi-particles and quasi-holes near the Fermi surface. 
In our case, the difference between quarks in different colors has been reflected by their 
propagators, we read the difference through calculating the quarks’ chiral condensate, but 
not trying to work out their masses. In the chiral limit, the chiral condensate disappears 
entirely in the color superconducting phase, it is not possible to investigate the influence of 
color breaking on quarks in different colors, so we will keep the current quark mass hnite in 
this paper. 

In the moderate baryon density region, people are interested in the question whether 
there exists a region where both the chiral symmetry and color symmetry are broken m 


r5| . In the chiral limit it was found that the existence of the mixed broken phase depends 


on the coupling constants Gs and Gum the quark anti-quark and diquark channels [|T^] |]T^ 


In the case of small ratio of Gd/Gs < 1 the calculations in the instanton model |T^ 


NJL model [0 and random matrix model [Q show a strong competition between the chiral 
and diquark condensates, i.e., where one condensate is nonzero the other vanishes. While 
for Gd/Gs > 1 the calculations in the random matrix model and NJL model show 
that there exists a region where both chiral (dynamical) symmetry and color symmetry are 
broken, and the chiral and diquark condensate coexist. The larger value of Gd/Gs is, the 
wider region of the mixed broken phase has been found in the random matrix model [l^ . 
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The presence of a small current quark mass will induce that chiral symmetry only restores 
partially and there will always exist a mall chiral condensate in the color superconducting 
phase, this phenomena had been called coexistence of chiral and diquark condensates in 
m- In order to differ from this coexistence of the diquark and the small chiral condensate 
induced by current quark mass, we will call the coexistence region of the diquark and large 
chiral condensate induced by large Gd/Gs before the chiral phase transition as mixed broken 
phase. 

In the coexistence region resulted by current quark mass, the chiral condensate is small 
comparing with the diquark condensate, and the role of the chiral condensate can be pre¬ 


scribed by the Anderson theroem [|T^, i.e., in this phase, the contribution of the chiral 
condensate to thermodynamic quantities becomes strongly suppressed, and one can calcu¬ 
late the diquark condensate neglecting the influence of the chiral condensate. 

In this paper, we will explain the existence of the mixed broken phase induced by larger 
Gd/Gs by analyzing the influence of the diquark condensate on the Fermi surface. In the 
mean-field approximation of the NJL model, the thermal system of the constituent quarks 
is a nearly ideal Fermi gas, and there is a sharp Fermi surface. The chiral symmetry begins 
to restore when the chemical potential is larger than the quark mass in the vacuum. When 
a diquark condensate is formed, the Cooper pair extends the Fermi surface, which induces 
the chiral symmetry restoring at a smaller chemical potential. The stronger the coupling 
constant in the diquark channel, the larger the diquark condensate and the smoother the 
Fermi surface. 

In the following, we briefly introduce the extended NJL model in Sec. II, then in Sec III, 
we evaluate the thermodynamic potential using the massive quark propogator. In Sec. IV, 
we get the gap equations and condensates. The numerical results and conclusions are given 
at the end. 
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II. THE EXTENDED NJL MODEL 


The choice of the NJL model |^] is motivated by the fact that this model displays the 
same symmetries as QCD and that it describes well the spontaneous breakdown of chiral 
symmetry in the vacuum and its restoration at high temperature and density. The model we 
used in this paper is an extended version of the two-flavor NJL model including interactions 
in the color singlet quark-antiquark channel and color anti-triplet diquark channel, which is 
not directly extended from the NJL model, but from the QCD Lagrangian |^- ||27|| . 

The importance of color 3 diquark degree of freedom is related to the fact that one can 
construct a color-singlet nucleon current based on it. Because the gluon exchange between 
two quarks in the color 3 channel is attractive, one can view a color singlet baryon as a 
quark-diquark bound state. And experimental data from pp collisions indicate the existence 
of this quark-diquark component in nucleons. 


The first attempt to investigate the diquark properties in the NJL model was taken in [28 


Starting from an NJL model for scalar, pseudoscalar, vector and axial-vector interactions of 
the {qq) x {qq) type and Fierz-transforming away the vector and axial-vector interactions, 
the scalar and pseudoscalar mesons, and diquarks can be obtained. However, this method 
could not get a consistent treatment of vector and axial-vector particles. 

The extended NJL model we used was derived directly from QCD Lagrangian 


Integrating out gluon degrees of freedom from QCD Lagrangian, and performing a local ap¬ 
proximation for the (nonperturbative) gluon propagator, one gets a contact current-current 
interaction. By using a special Fierz-rearrangement [^, one can completely decompose the 
two-quark-current interaction term into ’’attractive” color singlet (qq) and color antitriplet 
{qq) channels. In this way, a complete simultaneous description of scalar, pseudo-scalar, 
vector, and axial-vector mesons and diquarks is possible, thus the extended NJL model 
including {qq) x {qq) interactions is completed by a corresponding {qq) x {qq) interaction 
part. 

In our present work, we only consider scalar, pseudoscalar mesons and scalar diquark. 
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and using the following Lagrangian density 


n = q{iYd^ - mo)q + Gsligg)^ + ighbTgy] + GD[{i^ee^-f5g){igee''-f5g^')], 


( 1 ) 


where g'" = Gg^, (f = g^G are charge-conjugate spinors, G = ^ 7 ^ 7 ° is the charge conjuga¬ 
tion matrix (the superscript T denotes the transposition operation), mo is the current quark 
mass, the quark held g = gta with i = 1,2 and a = 1, 2, 3 is a havor doublet and color triplet, 
as well as a four-component Dirac spinor, r = (r^,r^,r^) are Pauli matrices in the havor 
space, where is antisymmetric, and are totally antisymmetric 

tensors in the havor and color spaces. 

In (|I|), Gs and Gd are independent ehective coupling constants in the scalar quark- 
antiquark and scalar diquark channel. The former is responsible for the meson excitations, 
and the later for the diquark excitations, which in principle can be determined by htting 
mesons’ and baryons’ properties in the vacuum. The attractive interaction in diherent 
channels in this Lagrangian will give rise to a very rich structure of the phase diagram. 
At zero temperature and density, the attractive interaction in the color singlet channel is 
responsible for the appearance of a quark anti-quark condensate and for the spontaneous 
breakdown of the chiral symmetry, and the interaction in the gg channel binds quarks into 
diquarks (and baryons), but is not strong enough to induce diquark condensation. As 
the density increases, Pauli blocking suppresses the gg interaction, while the attractive 
interaction in the color anti-triplet diquark channel will induce the quark-quark condensate 
around the Fermi surface which can be identified as a superconducting phase. 

After bosonization |^, one can obtain the linearized version of the model ([^) 


C = q{iYd^ - mo)g - g[a + * 7 ^ 7 ?)^ - ee^-i^g) - 


a2 + 7f2 

4 G 5 

with the bosonic helds 




AG 


( 2 ) 


D 


A® ~ ig‘"£e^75g. A*'’ ^ igee'^j^g^, c ^ gg, vf ~ ig7^rg. 


. *b 


C 


( 3 ) 
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Clearly, the a and tt fields are color singlets, and the diqnark fields and are 
color antitriplet and (isoscalar) singlet nnder the chiral SU{2)l x SU{2)r gronp. a ^ 0 and 
A^ 7 ^ 0 indicate that chiral symmetry and color symmetry are spontaneonsly broken. We 
assnme that the two condensates coexist, i.e., 

(T 7^ 0, vf = 0; 

a^ = a2 = o,aVo, (4) 

Here it has been regarded that only the first two colors participating in the condensate, 
while the third one does not. In the later expressions, we will simply nse A = A^. 

The real vacnnm will be determined by the minimnm of the thermodynamic potential at 
T = 0 and /i = 0, and the minimnm of the thermodynamic potential at any T, fi determines 
the stable state at that point. 


III. PARTITION FUNCTION AND THERMODYNAMIC POTENTIAL 

A. Nambu-Gorkov formalism 

The partition fnnction of the grand canonical ensemble can be evalnated by nsing stan¬ 
dard method IBDIl BTIl, 


Z = N' 


exp{ I dr J d^x {C + /xgqo*?)}, 


( 5 ) 


where /i is the chemical potential, and /? = 1/T is the inverse of temperatnre T. 

According to the mean field approximation Eq. (^), we can write the partition fnnction 
as a prodnct of three parts. 


7—7 7 7 

^ ^const-^qi^2'^Q3' 


( 6 ) 


The constant part is 


2const = N'exp{- dr j d^x [ 


cr^ A* A., 


( 7 ) 
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For the quarks in the hrst two colors (named as ’’the hrst two quarks” in the following) 
Q = ?i ,2 participating in the quark condensate, one has 

^qi ,2 = J[dQ][dQ]exp{ dr j d^x -m + n)Q + 

-m- ii)Q^ + ^QA"Q^ + ^Q^A+Q]}. (8) 

Here we have introduced the constituent quark mass 


and dehned A^ as 


m = mo + cr, 


( 9 ) 


A = —iAee^75, A’*' = —iA*£e%5 


( 10 ) 


with the relation A+ = 7 °(A“)l 7 °. 

For the quark in the third color (named as ’’the third quark” in the following), which is 
not involved in the diquark condensate, one has 


Zqz = j [dg 3 ][dg 3 ]exp{^ dr j d?x -m + -m- /i)gf]}. (11) 

Introducing the 8-component spinors for the third quark and the hrst two quarks, re¬ 
spectively 


^3 


(a ^ 

q3 

„c 

\Q 3 J 


^3 = (?3 qs), 


( 12 ) 


^ = 




, ^ = (Q Q^), 


and using the Fourier transformation in the momentum space. 




^-i{uJnT-p-x) 


""’qip), 


(13) 


(14) 


where V is the volume of the thermal system, we can re-write the partition function Eqs. 
(§) and (|TID in the momentum space as 







^ 91,2 ~ 


1 

exp{-E 

n,p 




(16) 


and 


^93 = /M^3]exp{^^ 

n^p 

= Det‘/2(/JG„-‘). (16) 

where the determinantal operation Det is to be carried ont over the Dirac, color, flavor 


and the momentum-frequency space. In Eqs. (^) and (W), we have defined the quark 
propagator in the normal phase 


Gn"' = 


/ 

GJ 

-1 

0 

\ 





-1 

V 

0 


Go- 

/ 


with 


( 17 ) 


[Col ‘ = (Po ± P)7o - 7 ■ P - ™ 


and po = iojn, and the quark propagator in the color-breaking phase 


f r^+1 ^ \ 


G"^ = 


V 


A+ 


Gn 


-1 


/ 


(18) 


(19) 


The Nambu-Gorkov propagator G(p) is determined from solving 1 = G ^ G, resulting in 


/ 


G = 


G+ 




G- 


( 20 ) 


with the components 


G^ = 


G 


1 -1 


- Gi]= , 


-1 


± _ 


■* = -GT A± Gj = -GJ A* G*. 


( 21 ) 


Here all components depend on the 4-momentum 
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In the case of the mass term m = 0, i.e., both the current quark mass tuq and the chiral 
condensate < a > are zero, the Nambu-Gorkov quark propogator has a simple form which 
could be derived from the energy projectors for massless particles [^] [^], while if there is 
a small mass term, the quark propagator can be expanded perturbatively around m = 0 , 
but its form is very complicated |^. In our case the quark mass term can not be treated 
perturbatively, we have to hnd a general way to deal with the massive quark propagator. 

Fortunately, we can evaluate a simple form for the massive quark propagator by using 
the energy projectors for massive particles. The energy projectors onto states of positive 
and negative energy for free massive particles are dehned as 


A,(p) = 1(1 ± 


( 22 ) 


where the quark energy Ep = y/fp + m^. Under the transformation of 70 and 75 , we can get 
another two energy projectors A±, 




which satisfy 


(23) 


7oA±(p)7i, = At(p), 75A±(p)75 = A±(p). 


(24) 


The normal quark propagator elements can be re-written as 


= 


7oA+ 


+ 


7oA 


Po + E^ Po-E^' 


with Ep = Ep± The propagator has four poles, i.e.. 


(25) 


Po = EEp , Po = tE^, 


(26) 


where the former two correspond to the excitation energies of particles and holes, and the 
latter two are for antiparticles and antiholes, respectively. 

The quark propagator including the diquark condensate’s contribution can be evaluated 


as 
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(27) 


G± = 


. Pi> - -Bp* 

P§ - Ef 


7o^+ + 


Po + E^ 
Pi - Ef 


70^-)(<^a/3 — ^a3<^/33)<5 




and 


, A± ~ A± ~ , 
pI - e% pI- ex 


(28) 


with EX‘^ = Ep"^ + A^. This propagator is very similar to the massless propagator derived 

in §. 

From the Nambu-Gorkov propagator, it is difhcnlt to obtain the mass of the qnark which 
participates in the diqnark condensate. The fonr poles of the Nambn-Gorkov propagator, 

i.e., 


Po = ±£^a, Po = =fBJ, 


(29) 


correspond to the excitation energies of qnasi-particles (qnasi-holes) and quasi-antiparticle 
(quasi-antiholes) in the color breaking phase. These quasiparticles are superpositions of 
particles and holes. 

We plot the excitation spectrum at p = SOOMeV as a function of Ep with different values 
of A//i in Fig. |^. A//i = 0 (the circles) correspond to the excitation spectrum in the 
normal phase. A/p = 0.2 (the squares), 0.5 (the stars) and 1 (the triangles) correspond 
to the excitation spectrum in the color superconducting phase, the black points are for the 
particles and the white points are for the holes, (a) is for the (quasi-)particles Ep{E~^) and 
(quasi-)holes —Ep{—E'^), and (6) is for the (quasi-)antiparticles —E^i—EX) and (quasi- 
)antiholes Ep{EX)- 

It can be easily seen from (a) that the quasi-particles and the quasi-holes mix particles 
and holes, which are called ’’Bogoliubons”. In the normal phase, to excite a pair of particle 
and hole on the Fermi surface does not need energy, while in the superconducting phase, 
to excite a pair of quasi-particle and quasi-hole at least needs the energy 2A when Ep = p. 
With increasing A, it is found that to excite a quasi-particle needs a larger energy, and the 
difference between the excitation energies sX Ep ^ p and that at Ep = p becomes smaller. 
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From ( 6 ), we see that to excite an antiparticle is mnch more difficnlt, and the diqnark 
condensate has little effects on the excitation spectrnm of (qnasi-)antiparticles and (qnasi- 
) antiholes. 


B. The calculation of InZgg 


For the third qnark which does not participate in the diqnark condensate, from Eq. (|T6|), 
we have 

lnZ„ = ll,i{Det(;3|G„]-‘) = ilii|Det(;3|G+]-')Det(;3[G„-]-')]. (30) 

Using the Dirac matrix, we first perform the determinant in the Dirac space. 


Det f3 [G^] ^ = Det (3 [{po + -j-p-m] 


/ 


= Det (3 


{po + p) -m a ■ p 


\ 


^ -3-p + 


= -;3" \(P« + P‘-E% 


(31) 


and in the similar way, we get 


Det (3 [Gq] ^ = -(3^ [(po - pf - E^]- (32) 

After performing the determinant in the Dirac space, we have 


Det f3 [G+] ^Det (3 [Gq] ^ = (3‘^[pl - {Ep + pf] (3‘^[pl - {Ep - pf], (33) 


Considering the determinant in the flavor, color, spin spaces and momentnm-freqnency 
space, we get the standard expression 

= N,Y. E{‘‘i(/3>o - (-Ep + pf]) + ln(/3>o “ (E, - p)=])}, (34) 

n f 

remembering that the color space for the third qnark is one-dimensional. 
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C. The calculation of InZ, 


91,2 


It is more complicated to evaluate the thermodynamic potential for the quarks partici¬ 
pating in the diquark condensate. From Eq. we have 


lnZqj 2 = -lnDet(/3G ). 

For a 2 X 2 matrix with elements A, B, C and D, we have the identity 


Det 


^ A 


C D 


= Det{-CB + CAC-^D) = Det(-5C' + DC-^AC). 


To prove the above equation, we have used 

/ 


A B 
C D 


^ ( {) B^ ^ 


C 0 


1 

B-^A 


C-^D ^ 


/ 


BC-^ AB 

DC-^ CB-^ 


^0 


B 0 


Replacing A, R, C and D with the corresponding elements of G we have 


Det(/3G-^) = /^^DetDi = f3^I)et[-A+A- + A+[G+]-^[A+]-i[Go ]-i] 

= yDetflj = yDet|-A+A- + [G„-]-‘[A-]-‘|GJl-‘A-]. 


(36) 


(36) 


(37) 


(38) 


Using the energy projectors A±, we can work out Di and D 2 as 


L>i = A^ + 75[7 o(Po - Ep )A+ + 7o(po + ^;^')A-]75[7o(Po - A;+)A+ + 7o(po + E^ )A_] 

= -U - Ef - A^)K_ + {pl - Ef - A^)Kl 

D 2 = -M - {E;f - A^)K + {pI - {E;f - A^)A_]. (39) 


Using the properties of the energy projectors, we can get 

D,D2 = U - {E;r - A^)] U - - A^)] = [pI - Ef][pl - Ef]. (40) 

With the above equations, Eq. (|3^) can be expressed as 

\nZ„, = lln[Det,3G-"] = A:An\fD^lj'-D 2 ] 

= l{lHti|^(2(pS - + TWn|^(^(p? - Bp)]} 

= 2/V/EE{1''1/3"(Po - -Ea')] + ln|/3"(P? “ Ef )|}. (41) 

n p 
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D. The thermodynamic potential 


The frequency summation of the free-energy 




( 42 ) 


can always be obtained by performing the frequency summation of the propagator 1 /(pg 


Ep). Differentiate Eq. (^) with respect to Ep\ 
d\nZf „ „ 1 


dE, 


= -2e,y: 


'r n Pd- Ej 


= I3[l - 2f(E,)], 


(43) 


where /(x) = l/(e^® + 1) is the usual Fermi-Dirac distribution function. Then integrating 
with respect to Ep, one can get the free-energy 


InZf = p[Ep + 2Tln(l + 


(44) 


With the help of the above expression, and replacing 

(Ep 




(27r)3’ 


we get the expressions 


InZga = 


(Ep 


^\E* + 2riii(l + e-”®'*) + E- + 2riii(l + )], 


(46) 


(46) 


= 


d^p 


[EX + 2Tln(l + + E^ + 2Tln(l + 


Finally, we obtain the familiar expression of the thermodynamic potential 


Q= -T 


\nZ 


a 


+ ^ ■ 2 ^^/ ( 1 ^ 1 '®!’ + + rin(l + ) 


d^p 


(47) 


1/ AGs ' ^Gd J (27r)3' 

+EX + 2Tln(l + e-^-^A) + + 2Tln(l + c-^-^a)]. 


(48) 
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IV. CONDENSATES AND GAP EQUATIONS 


A. Condensates 


With the Nambu-Gorkov quark propagator Eqs. ( 0 ) and (]28|), the diquark condensate 
is generally expressed as 

< >= {iT Y) [ (49) 

From general consideration, there should be eight scalar diquark condensates [0. In the 
case of the NJL type model, the diquark condensates related to momentum vanish, and there 
is only one independent 0+ diquark gap with Dirac structure F = 75 for massless quark, and 
there exists another O’*" diquark condensate with Dirac structure F = 7075 at nonzero quark 
mass. In our paper, we assume the contribution of the diquark condensate with F = 7075 is 
small, and only consider the diquark condensate with F = 75 . 

Performing the Matsubara frequency summation and taking the limit T —0, we get the 
diquark condensate at hnite chemical potential 


J- J- 

< (f>= —2ANcNf [ - —-1-17]. 

^ M 27i)^^ 2EX 2Et' 


(50) 


(27i; 

For the third quark, its chiral condensate can be evaluated by using the quark propagator 
in the normal phase, 

d^p 


< * 9 ’ >= 


(27r) 




(51) 


while for the quarks participating in the diquark condensate, the chiral condensate should 
be evaluated by using the quark propagator in the color breaking phase, 

< qiq^ >= -iT Y J 


We have no explicit mass expression for the hrst two quarks which participate in the 
diquark condensate, the influence of diquark condensate has been reflected in the quark 
propagator. The difference between the hrst two quarks which participate in the diquark 
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condensate and the third quark which does not participate in the diquark condensate can 
be read from their chiral condensates and can be dehned as 

S =< qiqi - < g 3<?3 (53) 

where 6 has the dimension of energy. In the case of chiral limit, the quark mass m decreases 
to zero in the color superconducting phase, and the influence of the diquark condensate on 
quarks in different colors vanishes. 

After performing the Matsubara frequency summation and taking the limit T —> 0, we 
get the expressions of the condensates at p 7 ^ 0 , 

< * 9 '’ > = imN, j - E,) - 1 |. 

/ d’^T) 1 

( 54 ) 

where 

1 

<= 2 ( 1 ^^) (55) 

are the occupation numbers for quasi-particles and quasi-antiparticles at T = 0. Correspond¬ 
ingly, 1 — Hp are the occupation numbers of quasi-holes and quasi-antiholes, respectively. 

We plot the occupation numbers for (quasi-)particles rip and (quasi-)holes I —rip in Fig. 
H (a), and the occupation numbers for (quasi-)antiparticles rip and (quasi-)antiholes 1 — n." 
in (b) as a function of Ep with respect to A//i = 0 (circles), 0.2 (squares), 0.5 (stars) and 1 
(triangles); the black and white points correspond to particles and holes, respectively. 

It is seen that the Fermi surface is very sharp in the normal phase A/p = 0, and becomes 
smooth when diquark condensate appears. The smearing is a consequence of the fact that 
the ’’Bogliubons” are superpositions of particle and hole states. The smearing of the Fermi 
surface induces the chiral symmetry restoring at a smaller chemical potential. The larger 
the diquark condensate is, the smoother the Fermi surface will be. 

From Fig. |^( 6 ), it can be seen that the occupation numbers for the (quasi-)antiparticles 
(antiholes) in the normal phase or color breaking phase are not sensitive to the magnitude 
of the diquark condensate. 
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B. Gap equations 


The two gap equations m and A can be derived by minimizing the thermodynamic 
potential Eq. (^) with respect to m and A, 

dVt dVt 


dm dA 


= 0 . 


(56) 


Taking into account the general expressions for the diquark condensate and chiral conden¬ 
sates, one can get the relations between the chiral gap m and the chiral condensate < qq >, 
i.e., 


m = mo -|- a, 


a = -IGs < qq >, 


(67) 


here the chiral condensate should perform summation in the color space 

< qq >= 2 < qiq^ > + < >; (58) 

and the relation between the diquark gap A and the diquark condensate < (f^^q >, i.e., 

A = -2Gd < >, (59) 


substituting Eq. (|50|) into the above equation, the gap equation for the diquark condensate 
Eq. (|59D in the limit of T —> 0 can be written as 

d?p 


1 = m^NfCD I 


(27r)3'2E 


1 1 


(60) 


V. NUMERICAL RESULTS 

In this section, through numerical calculations, we will investigate the phase structure 
along the chemical potential direction, analyze the competition mechanism between the 
chiral condensate and diquark condensate, and discuss the influence of the color breaking 
on the quarks in different colors. 
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Before the numerical calculations, we should hx the model parameters. The current 
quark mass mo = 5.5MeV, the Fermion momentum cut-off Aj = 0.637GeV, and the coupling 
constant in color singlet channel Gs = 5.32GeV~^ are determined by htting pion properties. 
The corresponding constituent quark mass in the vacuum is taken to be m(/i = 0) = 
330MeV. The coupling constant in the color anti-triplet channel Go can in principle be 
determined by htting the nucleon properties. In Go/Gs — 2.26/3 was chosen by htting 
the scalar diquark mass of ~ 600MeV to obtain a realistic baryon mass in the order of 
~ 900MeV. In our case, to investigate the inhuence of diquark condensate on the chiral 
phase transition, we will set Go/Gs = 0, 2/3,1,1.2,1.5, respectively. 


A. Phase structure at zero temperature 


First, we investigate the phase structure along the chemical potential direction with 
respect to diherent magnitude of Go/Gs- In the explicit chiral symmetry breaking case, we 
dehne the point at which the chiral condensate has maximum change as the critical chemical 
potential for the chiral phase transition, and the point at which the diquark condensate 
starts to appear as the critical chemical potential pa for the color superconductivity phase 
transition. 

The two gaps m (white points) and A (black points) determined by Eqs. (|T^) and (|5^ 


are plotted in Fig. ^ as functions of /i with respect to different Go/Gs = 2/3,1,1.2,1.5 
in (a), (5), (c) and (d), respectively. In Fig. they are plotted as functions of the scaled 
baryon density Ub/no, where uq is the normal nuclear matter density. 

In Fig. 1^ (a), i.e., in the case of Go/Gs = 2/3, we see that in the region where the 
constituent quark mass keeps its value in the vacuum m(/i = 0), the diquark condensate 
keeps zero. The chiral phase transition and the color superconductivity phase transition 
nearly occur at the same chemical potential ~ /ia = 340MeV. The two phase transitions 
are of hrst order. In the explicit chiral symmetry breaking case, there is a small chiral 
condensate in the color superconductivity phase p > This phenomena has been called 
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the coexistence of chiral and diquark condensate in |^. In this coexistence region, the 


chiral condensate is small and can be described by the Andersom theorem 


In (6), with Gd/Gs = 1, the diquark condensate starts to appear at /xa = 298MeV, then 
chiral symmetry restores at = 304.8MeV. Both /ta and are smaller than those in the 
case of Gd/Gs = 2/3. In the region from /ta to both chiral and color symmetries are 
broken, which is called the mixed broken phase. The diquark gap increases continuously 
from zero to 82MeV in this mixed broken phase, and jumps up to 152MeV at the critical 
point /i^. The chiral symmetry phase transition is still of the hrst order, and the jump of 
the diquark gap can be regarded as the influence of the hrst order chiral phase transition. 

In (c) and (d), we see that with increasing of Gd/Gs, the diquark condensate starts to 
appear at smaller /ia, and chiral symmetry restores at smaller while the width of the 
region of mixed broken phase, becomes larger. This phenomena has also been found 

in Ref. |]T^ in the random matrix model. It is also seen that with increasing of Gd/Gs, the 
hrst order phase transition of chiral symmetry restoration becomes the second order one, 
where the critical point /i^ is still dehned by the maximum change of quark mass m. 

We summarize the phase structure along the chemical potential direction: 1) when /i < 
/xa; chiral symmetry is broken; 2) in the region from /xa to both chiral and color 
symmetries are broken, and 3) when /i > /x^, chiral symmetry restores partially and the 
phase is dominated by color superconductivity. The phase transition from chiral broken 
phase to the mixed broken phase is of second order, and the phase transition from the 
mixed broken phase to color superconductivity phase is of hrst order for small Gd/Gs, and 
of second order for large Gd/Gs- The phase structure depends on the magnitude of Gd/Gs- 
If Gd/Gs is very small, the diquark condensate will never appear; If Gd/Gs < 1 but not 
too small, there will be no mixed broken phase, the chiral phase transition and the color 
superconductivity phase transition occur at the same critical point = /xa, and the two 


phase transitions are of hrst order like in |10 
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B. The Competition between the chiral and diquark condensate 


With increasing of Gd/Gs, the diquark condensate starts to appear at smaller and 
chiral symmetry restores at smaller /i^, while the width of the mixed broken phase, —//a, 
increases. 

In order to understand the competition mechanism and explicitly show how the diquark 
condensate influences the chiral phase transition, we plot the constituent quark mass m and 
the diquark gap A as functions of /i for different values of Gd/Gs in Fig. 

In the case of Gd/Gs = 0, only chiral phase transition occurs, the thermal system in the 
mean-field approximation is nearly a free Fermi gas made of constituent quarks. In the limit 
of T = 0, there is a very sharp Fermi surface of the constituent quark. When the chemical 
potential is larger than the constituent quark mass in the vacuum, the chiral symmetry 
restores, and the system of constituent quarks becomes a system of current quarks. When 
a diquark gap A forms in the case of Gd/Gs ^ 0, smooth the sharp Fermi surface of 

the constituent quark. In other words, the diquark pair lowers the sharp Fermi surface, and 
induces a smaller critical chemical potential of chiral restoration. 

In TABLE I, we list the chemical potentials /ia, at which the diquark gap starts ap¬ 
pearing, and /i^, at which the chiral symmetry restores, for different values of Gd/Gs- 

= 345.3MeV is the critical chemical potential in the case of Gd/Gs = 0. A^ is the value 
of diquark gap at jj,^, if there is a jump, it is the lower value. 

We see that for larger Gd/Gs, the diquark condensate appears at a smaller chemical 
potential pA, and the chiral phase transition occurs at a smaller critical chemical potential 
/i^, the gap of the diquark condensate A^ at becomes larger, and the region of mixed 
broken phase becomes wider. 

We assume the relation between A^ and as 

hx = Ff - (61) 

and the relation between A^ and /ia as 
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Aa = fiF - yAj. 


( 62 ) 


In TABEL I, we listed the values of x = {fip — fiy.)/A^ and y = {yPp — for 

different Gp/Gs- It is found that x is almost Gp/Gs independent and equal to 1 / 2 . As for 
y, it is larger than 1/2 and increases with increasing of Gd/Gs- From the Eqs. (^) and 
we have the relation 

hx-hA = ( 2 /-1/2)A^ (63) 

for the mixed broken phase. With increasing of Gd/Gs, y and increase, and then the 
width of the mixed broken phase, ~ hA, becomes larger. 

Now we turn to study how the chiral gap influences the color superconductivity phase 
transition. Firstly, we change the constituent quark mass in the vacuum from 330MeV to 
486MeV. To £t the pion properties, the coupling constant in the quark-antiquark channel 
is correspondingly increased from Gs to I. 2 G 5 . We plot the diquark gap as a function of y 
in Fig. I^a for the two vacuum masses and for Gd/Gs = 2/3,1,1.2,1.5. We hnd that for the 
same Gd/Gs, the diquark gap starts to appear at much larger chemical potential /xa when 
the vacuum mass increases from 330MeV to 486MeV. 

Then we withdraw the quark mass, i.e., taking m = 0 even in the vacuum. We plot 
the diquark gap as a function of fi in Fig. ^ h for m(/i = 0) = 0 and 330MeV and for 
Gd/Gs = 2/3,1,1.2,1.5. We hnd that for any Gd/Gs, the diquark condensate starts 
to appear at a much smaller chemical potential /ia for m(/i = 0 ) = 0 compared with 
m(/i = 0) = 330MeV. 

From Fig. a and b, we can see that quark’s vacuum mass only changes the critical point 
of color superconductivity /ia, the diquark gaps for different quark’s vacuum mass coincide in 
the overlap region of color superconductivity phase, where chiral symmetry restores partially. 

From the inhucence of the diquark gap on the chiral phase transition and the inhuence 
of the chiral gap on the color superconductivity phase transition, it is found that there does 
exist a strong competition between the two phases. The competition starts at pa and ends 
at /i^. We call the mixed broken phase as the competition region, which becomes wider with 
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increasing of Gd/Gs- This competition region is the result of the diquark gap smoothing 
the sharp Fermi surface of the constituent quark. If the attractive interaction in the diquark 
channel is too small, there will be no diquark pairs, the system will be in the chiral breaking 
phase before /i^, and in the chiral symmetry restoration phase after If the attractive 
interaction in the diquark channel is strong enough, diquark pairs can be formed and smooth 
the sharp Fermi surface, and induce a smaller critical chemical potential of chiral symmetry 
restoration. 


C. The influence of color breaking on quarks’ properties 

Finally, we study how the diquark condensate influences the quark properties. 

In the normal phase, the quarks in different colors are degenerate. However, in the color 
breaking phase, the hrst two quarks are involved in the diquark condensate, while the third 
one is not. 

The quark mass m appeared in the formulae of this paper is the mass for the third 
quark which does not participate in the diquark condensate. We have seen from Fig. ^ (a) 
that the diquark condensate influences much the quark mass m in the competition region 
Pa < /i < Px- color breaking phase, i.e., when fi > pA, the quark mass m in 

different cases of Go decreases slowly with increasing p, and reaches the same value at 
about p = 500MeV. 

The difference of the chiral condensates for quarks in different colors 6 dehned in Eq. 
(^) is shown in Fig. ^ as a function of the chemical potential fi with respect to Gd/Gs = 
2/3,1,1.2,1.5. It is found that in any case S is zero before pa, then begins to increase at 
Pa and reaches its maximum at Px, and starts to decrease after p > p^, and approaches to 
zero at about p = 500MeV, when p > 500MeV, S becomes negative. With increasing Gd, 
S{fi^) increases from IMeV for Gd/Gs = 2/3 to 13MeV for Gd/Gs = 1-5. Comparing with 
the magnitude of the diquark condensate, 5 is relatively small in the color superconductivity 
phase. 
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VI. CONCLUSIONS 


In summary, in an extended NJL model and considering only the attractive interactions 
in the O’*" color singlet quark anti-quark channel and color anti-triplet diquark channel, 
the Nambu-Gorkov form of the quark propogator has been evaluated with a dynamical 
quark mass. The Nambu-Gorkov massive propagator makes it possible to extend the Green- 
function approach to the moderate baryon density region, and the familiar expression of the 
thermodynamic potential has been re-evaluated by using the massive quark propagator. In 
this paper, we have neglected another scalar diquark condensate with Dirac structure 7075 , 
which is assumed to be small. 

The phase structure along the chemical potential direction has been investigated. The 
system is in the chiral breaking phasei before pa, in the color superconducting phase after 
and the two phases compete with each other in the mixed broken phase with width — Pa- 
The width depends on the magnitude of Gd/Gs- If G^/Gs is small, the width is zero, and 
the chiral phase transition and the color superconductivity phase transition occur at the 
same chemical potential. The width increases with increasing of Gd/Gs- The competition 
mechanism has been analyzed by investigating the influence of the diquark condensate on 
the sharp Fermi surface. The diquark condensate smoothes the sharp Fermi surface, and 
induces the chiral phase transition occuring at a smaller chemical potential, and the diquark 
can be formed more easily. The phase transition is of second order from the chiral symmetry 
broken phase to the mixed broken phase, and the phase transition from the mixed broken 
phase to the partial chiral symmetry phase, i.e., the color superconductivity phase is of first 
order for small Gd/Gs, and is of second order for large Gd/Gs- 

The influence of the diquark condensate on the properties of quarks in different colors 
has also been investigated. It is found that the difference of the chiral condensates between 
quarks in different colors induced by the diquark condensate is very small. 
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1 



FIG. 1. The excitation spectrum at ^ = 500MeV for (quasi-)particles and (quasi-)holes in (a) 
and for (quasi-)antiparticles and (quasi-)antiholes in (b) as a function of Ep with different values 
of A/fj,, A/fi = 0 (circles), 0.2 (squares), 0.5 (stars) and 1 (triangles). The black and white points 
are for the particles and holes, respectively. 
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FIG. 2. The occupation numbers for (quasi-)particles and (quasi-)holes 1 — n+ in (a), and 
for (quasi-)antiparticles n~ and (quasi-)antiholes 1 — n“ in (b) as a function of Ep with respect 
to A/fi = 0 (circles), 0.2 (squares), 0.5 (stars) and 1 (triangles). The black and white points 
correspond to particles and holes, respectively. 
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FIG. 3. The two gaps m (white points) and 


fi for Gd/Gs = 0, 2/3,1,1.2,1.5, respectively. 
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FIG. 4. The two gaps m (white points) 
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FIG. 5. The gaps m in (a) and A in (6) as a function of the chemical potential /x with respect 


to Gd/Gs = 0, 2/3,1,1.2,1.5, respectively. 
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FIG. 6. The influcence of chiral gap on the color superconductivity phase transition in the 


case 


of miji = 0) = 486MeV in (a) and m(/U = 0) = 0 in ib). 
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FIG. 7. The difference of the chiral condensates for quarks in different colors d as a function of 
the chemical potential // with respect to Gd/Gs = 0, 2/3,1,1.2,1.5, respectively. 
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Gd/Gs 

/iA(MeV) 

/^x(MeV) 

^J^x - liA(MeV) 

Ax(MeV) 

ifJ-F - f^x)/^X 

(/^F - Ma)/Ax 

1 

298 

304.8 

6 

82 

0.49 

0.57 

1.2 

242 

266 

24 

162 

0.49 

0.64 

1.5 

70 

190 

120 

310 

0.50 

0.89 


TABLE I. The Gd dependence of chemical potentials and iJp = 345.3MeV. 
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